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Abstract 

We propose a free falling particle in an AdS space as a holographic model of local quench. 
Local quenches are triggered by local excitations in a given quantum system. We calculate the 
time-evolution of holographic entanglement entropy. We confirm a logarithmic time-evolution, 
which is known to be typical in two dimensional local quenches. To study the structure of 
quantum entanglement in general quantum systems, we introduce a new quantity which we 
call entanglement density and apply this analysis to quantum quenches. We show that this 
quantity is directly related to the energy density in a small size limit. Moreover, we find a 
simple relationship between the amount of quantum information possessed by a massive object 
and its total energy based on the AdS/CFT. 



Contents 

1 Introduction Q 

2 Holographic Local Quenches as Falling Particles 

2.1 A Falling Massive Particle in AdS 

2.2 Einstein Equation 

2.3 Coordinate Transformation from Global AdS to Poincare AdS 

2.4 Back-reacted Metric for a Falling Massive Particle 

2.5 Holographic Energy Stress Tensor 

3 Perturbative Analysis of Holographic Entanglement Entropy under Local Quenches 

3.1 Perturbative Calculations of HEE under Local Quenches 

3.2 Explicit Calculations of HEE under Local Quenches 

3.3 Small Subsystem Limit: An Analogue of the first law of thermodynamics 

3.4 Large Subsystem Limit 

4 Exact Holographic Entanglement Entropy for 2d Local Quenches 

4.1 Symmetric Intervals 

4.1.1 Case 1: M < 

4.1.2 Case 2: M > B? 

4.1.3 Final Results 

4.1.4 Thermal Local Quenches 

4.2 General Formulation 

4.3 An Interval writh An Excited End Point 

4.3.1 Semi Infinite Limit i — > oo 

4.4 HEE for General Intervals 

5 Entanglement Density 

5.1 Definition of Entanglement Density 

5.2 Strong Subadditivity and Positivity of Entanglement Density 

5.3 Conservation Law 

5.4 Entanglement Density for Global Quenches 

5.5 Entanglement Density for Holographic Local Quenches 

6 Quantum Information, Thermodynamics and Gravity 

6.1 Total Amount of Quantum Information 

6.2 Entanglement Renormalization and Origin of Gravitational Force 

7 Conclusions and Discussions 

A Perturbative Construction of Gravity Duals 

A.l Perturbation Theory 

A. 2 General Solutions to Einstein Equation 

A. 3 Special Solutions: a— ^-O 



1 



1 Introduction 



Quantum entanglement provides us with a powerful method of investigating various quantum states 
and classifying their quantum structures. Moreover, quantum entanglement is a useful tool when 
we would like to study excited quantum systems which are far from thermal equilibrium, for ex- 
ample, those under thermalization processes. Even though in such systems we cannot define the 
thermodynamical entropy and temperature etc, we can always define the entanglement entropy 
(refer to e.g. the reviews [H [21 [3l HI [5]). 

Quantum quenches provide us with ideal setups to study thermalizations of quantum systems, 
which can be realized even in real experiments such as cold atoms. They are unitary evolutions 
of pure states triggered by sudden change of parameters such as mass gaps or coupling constants 
[U [71 [8]. By using the AdS/CFT O [lOl [11], we can relate this problem to the dynamics of 
gravitational theories [12l [131 [HI [E] • 

One typical class of quantum quenches is called global quenches and they occur from homo- 
geneous changes of parameters [6|. Holographic duals of global quenches are dual to black hole 
formations as have been discussed in |121I14|. A useful quantity to understand how thermalizations 
occur is the entanglement entropy ^ . The holographic entanglement entropy (HEE) [161 llTl [5] 
has been calculated for various quantum quenches [131 \T8\ [19] . Refer also to [20] for a computation 
of HEE for a stationary system which is described by an excited pure state dual to an AdS plane 
wave. In general, non-local probes such as the entanglement entropy are useful to measure the 
thermalization time. On the other hand, local quantities such as an expectation value of energy 
momentum tensor shows an immediate thermalization and is not suitable to see if a given system is 
completely thermalized [18]. In this sense, the entanglement entropy can serve as a non-equilibrium 
substitute of thermodynamical entropy. 

Another class of quantum quenches is the local quench and this is triggered by a shift of 
parameters within a localized region or simply at a point. One of the aims of the present paper 
is to provide a simple construction of holographic dual for local quench and is to calculate the 
holographic entanglement entropy. A local quench shows how localized excitations in a given 
quantum system propagate to other spatial regions. Local quenches have been studied in two 
dimensional CFTs [3 [8]. However, local quenches in higher dimensions have not been understood 
well. This partially motivates us to study the local quenches in AdS/CFT, which often allows us 
higher dimensional calculations. 

We will argue that a simple holographic description of a system just after the local quench is 
a free falling particle-like object in an AdS space. It is pulled into the horizon of AdS space due 
to the gravitational force and this is the reason why we observe the non-trivial time-dependence 
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of entanglement entropy. Therefore this problem is deeply connected to a fundamental question: 
what is the non-gravitational (or CFT) counterpart of gravitational force via holography ? We will 
suggest an intuitive answer to this question in the end. 

The time evolution of quantum entanglement under local quenches is more complicated than 
the global quenches because it is inhomogeneous. To understand its structure clearly we introduce 
a new quantity which we call entanglement density. It is defined by taking the derivatives of the 
entanglement entropy Sa with respect to the positions of two boundary points of the subsystem A. 
This quantity counts the number of entangled pairs at a given position. The strong subadditivity 
guarantees that this quantity is always positive. As we will see this analysis reveals the detailed 
structure of quantum entanglement under local quenches as well as global quenches. 

One more motivation to study local quenches is to estimate the amount of quantum information 
possessed by a massive object or radiations. We will employ the entanglement entropy for local 
quenches in order to measure the amount of information included in a localized excited lump. We 
will evaluate this quantity by using the holographic entanglement entropy (HEE) and obtain the 
simple conclusion that it is given by the total energy of the object times its size up to a numerical 
factor. 

The paper is organized as follows: In section two, we explain our holographic setup of local 
quench using a free falling particle in AdS^+i. We calculate the holographic energy stress tensors 
in this model. In section three, we compute the holographic entanglement entropy for d = 2,3,4 
assuming that the back-reaction due to the falling particle is very small. In section four, we perform 
an exact analysis of holographic entanglement entropy for d = 2. In section 5, we introduce a 
new quantity which we call entanglement density and we investigate the evolution of quantum 
entanglement structures under local quenches by using this. In section 6, we study the relation 
between the amount of information of an object and its total energy using the results in the 
previous sections. We also interpret our results of local quenches using the idea of the entanglement 
renormalization and discuss the holographic interpretation of gravitational force. In section 7 we 
summarize our conclusions and discuss future problems. In appendix A, we show an explicit 
perturbative calculation of back-reactions due to the falling particle. 

2 Holographic Local Quenches as Falling Particles 

Local quenches in quantum systems are triggered by a sudden local change of the Hamiltonian at 
a specific time. One typical class of examples will be joining two separated semi-infinite systems 
at each endpoint as studied in [71 [8] (see the upper picture in Fig{T]). When this quench process 
happens, an interaction between two endpoints is instantaneously introduced. From the viewpoint 
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Figure 1: Setups of local quenches. The upper picture describes a process of jointing two systems 
which are defined on semi-infinite lines. The lower one describes localized excitations on an infinitely 
extended system. We define the parameter a which measures the size of excited region at the 
beginning of the local quench. 

of the new Hamiltonian, an locally excited state is generated just after this local quench. Therefore 
we can generally characterize a local quench by local excitations. These excitations will propagate 
to other regions under the time-evolution. 

Now we would like to construct gravity duals for local quenches via the AdS/CFT. Even though 
we have not found a simple gravity dual of the original model i.e. joining two CFTs, it is not 
difficult to find a holographic model for local excitations. It is given by a falling massive particle 
in a Poincare AdS space (see Figl2]). At t = 0, the particle is situated near the AdS boundary and 
its back-reaction to the metric is highly localized near the particle. Under its time evolution it falls 
into the AdS horizon and its back-reaction spreads out. In the dual CFT, at t = 0, the excitations 
are concentrated in a small localized region (see the lower picture in Fig{T]), whose radius is defined 
to be a. Therefore we can regard this state at t = as the one just after the local quench. Later 
the excitations expand at the speed of light as we will see e.g. from its holographic energy stress 
tensor later. In this way, we can regard this setup as a gravity dual of local quench. 

We will employ the beautiful construction of back-reacted solutions found by Horowitz and 
Itzhaki in [JJ . The basic idea is to start from a black hole in a global AdS space and map it into a 
Poincare AdS by the coordinate transformation. This leads to a falling black hole solution. Though 
we are thinking of a massive particle with a finite size or equally a star, instead of a black hole, the 
asymptotic solution which is outside of the star is the same as that for a black hole as usual. 
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Figure 2: A falling massive particle in AdS and the calculation of holographic entanglement entropy 
for two different choices of the subsystem A. It is clear from this picture that the back reaction 
due to the falling particle gets significant when / = z(t) in the left picture and t = in the right 
one because the particle is on top of ^a- 



2.1 A Falling Massive Particle in AdS 

Consider a d + 1 dimensional AdS space (AdS(i_|-i) in the Poincare coordinate 



(2.1) 



The radius of AdS is defined to be R and the coordinate of AdS is represented by = {z, t,xi,- ■ 

■,Xd-.l). 

In this AdS space, we introduce a massive object (mass m) with a very small size which is 
larger than the Schwartzschild radius. Its motion in the AdS space is described by the trajectory 

= X^^ir). In general, the action of a particle with mass m in a spacetime defined by the metric 
g^,i, is given by 



-m 



dr / dx'^+i5(rf+i)(x/^ _ X^'{T))^-g^^{x) ■ ^rX^'{T) ■ drX'^ir). 



(2.2) 



We assume that the particle is situated at X* = and we gauge fix by setting X^{t) = r. Then 
the trajectory is specified by the function X^[t) = z[t). In the pure AdS background (12. ip . the 
action looks like 

^ ^' (2.3) 



-mR I dt 



z{t) ■ 

The solution to the equation of motion derived from ()2.3p is given by 



Z(t) = x/(t-to)2 + 



(2.4) 
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where 'ind a are integration constants. Below we will set = by using the time translation 
invariance. When t < the particle moves from the horizon to the boundary. It reaches z = a at 
t = 0. Later {t > 0), it again falls into the horizon as depicted in Figl2j Thus the energy of the 
particle in the AdS space is calculated as 

^ mi? 

E= . (2.5 

a 

2.2 Einstein Equation 

The gravity action coupled to the massive particle reads 

Stot = j dx''+^^{n - 2A) + S-p , (2.6) 

where the cosmological constant is given by A = — -^r^ and Gn is the Newton constant. 
The equation of motion becomes 

T^^" - ^g^"'^ + ^g'"' = T^", (2.7) 
where the bulk energy-stress tensor T^'^ is given by 

We will show a direct perturbative calculations of this back-reaction in appendix A. However, 
below we will take a different step in order to analytically construct the back- reacted solutions. See 
the paper |22] for analytical calculations of back-reactions to a scalar field in an AdS space. Refer 
also to [23] for a more extensive analysis and a relation to expanding qluon plasmas, where the 
back-reacted solutions are called conformal solitons (see [23] for spacetime structures of conformal 
solitons) . 

2.3 Coordinate Transformation from Global AdS to Poincare AdS 

Now consider the global AdS^+i space defined by the metric 

ds^ = -(i?2 + r^^dr' + + r^dnj_,. (2.9) 

_|_ 



We can show this is (locally) equivalent to the Poincare AdS^+i space ()2.ip via the following 
coordinate transformation: 

R2e^ ^e-^(z^ + -t'') 



+ cos r 



1^2 , 



2z 

Rt 



\/ R"^ + r"^ sin r = — , 

z 

rni = ^ (i = l,2,---,d-l), 

-W + e-/^(z^ + x^-t^) 
rild = . (2.10) 
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Here, the coordinate of S'^~^ is described by (Oi,il2,- • -j^d) such that Yli=ii^i)'^ — 1- Also 
we defined x'^ = ^2^=1^1- The arbitrary constant /3 is introduced for the later purpose, which 
corresponds to the boost transformation of S0{2,d) symmetry. If we set /3 = 0, ()2.10p is reduced 
to the standard one which can be found in e.g. 

2.4 Back-reacted Metric for a Falling Massive Particle 

In the global coordinate, we can consider a static particle situated at r = 0. Following the idea in 
|21] . we would like to map it into the Poincare AdS. After the coordinate transformation (|2.10p . 
its trajectory is mapped into 

Xi = 0, z'^ -t^ = R^e'^^. (2.11) 
Thus this corresponds to the previous trajectory (j2.4p with the identification 

a = Re^. (2.12) 

The back-reacted geometry outside of the massive object is obtained from the AdS black hole 
solution [25]: 



ds' = - r' + R'--iz^]dT'+ , + r'dnti- (2.13) 



M \ 2 R^dP 

Note that in the AdSa case (d = 2), the solution (j2.13p for M < B? \s not a black hole solution but 
a solution with a deficit angled] The mass parameter M in ()2.13p is related to the mass m of the 
particle via 

8r(d/2) GnR^' ^ ' 

Therefore, we can find the back-reacted metric by performing the coordinate transformation 

(12.10p to the metric (I2.13p . This can be done in a straightforward manner by noting 

1 



d 

dT^ = dicosrf + d{smTf, dnl_^ = J^i^^if- (2.15) 
2.5 Holographic Energy Stress Tensor 

One way to understand the time evolution of the CFT^ state dual to the falling particle in AdS^+i, 
is to calculate the holographic energy stress tensor. For this purpose, it is useful to employ the 
Fefferman-Graham gauge of the coordinates given by the expression 

d.2 = r2 . dz^ + g^,ix z)dx^dx' ^ 
z^ 



^We do not have to worry about the singularity because we replace the region near it with a star solution. 



7 



where x"" = {t,Xi). We are considering the case where the boundary metric gab{x,0) coincides with 
the flat Minkowski metric rjab as in the Poincare AdS. Then near the AdS boundary 2 = 0, gab{x, z) 
behaves hke 

gab{x, Z) = r]ab + tabZ'' + 0{z''+^). (2.17) 

In this setup, the holographic energy stress tensor [26] is calculated from the formula: 

Note that the metric we find from the coordinate transformation (|2.10p of (|2.13p is not in the 
form of the Fefferman-Graham gauge (j2.16p . Thus we need to perform a coordinate transformation 
further to achieve this gauge so that we can employ (j2.18p . It is also useful to define the light-cone 
coordinate u = t — p and v = t + p, where p = xi for d = 2 and p = \J X] ti -^"i d > 2. Finally 
we obtain the following energy stress tensor Tab for d = 2, 3,4: 

d = 2 (AdSs) : 

87rG^i?(n2 + a2)2' ^™ " 87rG^i?(t;2 + a2)2 ' ^™ " 0- (2-19) 



d = 3 (AdS4) : 



3Ma^ „ Ma^ 



= - Tuv — 

J-vv — 5 — 1 J^ee — 3 3"i vz.zuj 

87rG7vi?(^^2 + a2)iVn2 + a2 87rGjv«(«2 + q2)I(^;2 + ^^2)1 

d = A (AdSs) : 

T — T 

J- 71.11. _ _ _ - i-*-? 



T 



■KGNR{a'^ + v?f{a'^ + v'^y ™ 2^Givi? (^2 + n2)2 (a2 + f;2)2 ' 



2^3' 



^ _ Ma4(ti-u)2 _ MaSin^ (9(n - f;)^ 

~ 47rGivi? + u^f {a"^ + v^f ' ^'^^ ~ AttGnR (a^ + u^f {a^ + i;2)2 ' (2-21) 

The angular coordinates 9 for d = 3 and (0, (p) for d = 4 are those of the polar coordinates of 
R^ : dp^ + p^dO"^ and i?'^ : + p'^{d9'^ + sin^ OdcjP'), respectively. For ti = 4, the above expression 
agrees with the one in |21j . 

Note that we can easily confirm the traceless condition 

TabTi"'' = 0, (2.22) 
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and the conservation law 

daT^'' = 0. (2.23) 
In particular, the energy density Tu in each dimension is given by 



d = 2: Tu 
d = 3: Tti 



Ma3 (p2 + + a2)2 + 2pH^ 



)2 

. 4 T 3(p^ + t^ + a^)^+Vt^ 

** VrGTvi? ((/92-t2_c,2)2 + 4c,2^2)3- ^^'^^^ 

By using these expressions we can confirm 

'd^-^xTu = — = E, (2.25) 
a 

for any d, by using the relation (j2.14p . Therefore the total energy agrees with that of the particle 
()2.5p as expected. 

To understand how the excitations propagate, it is useful to look at the time evolution of Tu- 
This is plotted in Figj3j We can observe a peak on the light-cone = as can be understood as 
the light-like propagation (or shock waves as in [2Tj) of the initial excitations at t = p = 0. 

The parameter a parameterizes the width of this peak and therefore it measures the size of 
the localized excitations. Notice that for d = 2, the heights of the two peaks are equal and stay 
constant due to the total energy conservation. We sketched the essence of this behavior in Figj4] 
for d = 2. In the zero width limit a — t- we find Tu gets delta-functionally localized: 



d = 


2 : 


Tu- 


^^{Kt + p) + Kt-p)), 




d = 


3 : 


Tu- 






d = 


4 : 


Tu- 


> + P^ + - ■ 


(2.26) 



3 Perturbative Analysis of Holographic Entanglement Entropy 
under Local Quenches 

The entanglement entropy Sa is defined as the von-Neumann entropy when we trace out the 
subsystem B, which is the complement of A. The subsystem A is an arbitrary chosen space- 
like region on a time slice. Therefore, in a time-dependent background, Sa_ depends on the time 
even if we fix the shape of the region A. This is expected to be an important quantity which 
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Figure 3: The profiles of the energy density Tu ior d = 2 (left), d = 3 (middle) and d = 4 (right) 
as a function of p and t within the range —5 < t < 5. We set a = R = M = Gn = 1- 



Ttt 




Entangled Pair 



Figure 4: A sketch of time evolution of energy density and entangled pairs for d = 2. The under- 
standing of detailed structure of quantum entanglement and entangled pairs is the main subject of 
this paper as will be studied in later sections. 



characterizes various non-equilibrium processes in quantum many-body systems such as quantum 
quenches [H El [7] . 

In holographic setups, we can calculate Sa by using the holographic entanglement entropy 
(HEE). In d+1 dimensional static gravity backgrounds, we can calculate Sa by the formula |16j 

Area(7A) 



Sa 



(3.1) 



where ^a is the d — 1 dimensional minimal area surface on a time slice which ends on the boundary 
of A i.e. dA. For spherically symmetric subsystems, we can manifestly prove (j3.ip in pure AdS 
spaces via the bulk to boundary relation [lOj as shown in [27]. In time-dependent backgrounds, 
there is no natural time-slice in the bulk AdS and we need to employ the covariant version of HEE 
|17j . This is given by redefining 7^ to be an extremal surface in the Lorentzian spacetime. 

Below we would like to compute the holographic entanglement entropy in the gravity dual of 
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the local quench obtained from the coordinate transformation of the back hole solution (|2.13|) as 
explained in the previous section. In this section we perform a perturbative calculation assuming 
that the back-reaction from the particle is very small, keeping only the leading term proportional 
to M. This allows us a relatively simple calculation which is applicable to any dimension d. 

3.1 Perturbative Calculations of HEE under Local Quenches 

Consider the first order perturbation of the metric 

5/.. = i°j+5?J+0(M2), (3.2) 

where g^^J represents the metric of pure AdS^+i ()2.ip : gl:}u is the leading perturbation due to 
the back-reaction, which is of order M in our case. We obtained g^j^J from the direct calculation 
explained in the previous section, though we will not write its complicated expression explicitly. 

If we know the extremal surface 7_a in the pure AdS^+i, then the perturbed area of an extremal 
surface can be found as 

AArea =^J d'^-^^VG(^Tr[G^^\G^^'^)-^]. (3.3) 

Here G is the induced metric on the surface 7^: 

(0) ^ dJO^dJC^jo) ^(1) ^ dX^dJT m 

where ^" (q = 1, 2, • • •, (i — 1) is the coordinate of the codimension two surface ja- The embedding 
function X^{^) is that of 'ja in the pure AdS^^+i. Notice the useful fact that we do not need to 
know the precise shape of ja in the perturbed metric to calculate 



3.2 Explicit Calculations of HEE under Local Quenches 

Consider the d = 3 case (AdS4) first. We take the subsystem A to be a round disk with the radius 
/, defined by + < P. The corresponding extremal surface (or equally minimal surface) is given 
by the half sphere [16] 



'P-xl-xl (3.5) 
We can take (^"^,C^) = (xi,X2) and then we find that the area density is given by 



iv^TY[G«(GW)-^] = ^, (3.6) 

2 /(e2/3i^4 + g-2/3(|2_i2)2 + 2i?2(2p2+i2_^2))3/2 



where p= \/T.t=i^'i- 
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The area perturbation is found by integrating (j3.6|) over x and y. By plugging this into (j3.ip 
we obtain the increased amount of HEE, denoted by A5a, compared with Sa for the ground state 
dual to the pure AdS: 

4G7vi?a/ y + - i^) + (a2 + ^2)2 i y' ^ ^ 

where we used the relation (j2.12|) . Note that the area law divergence [25] is canceled in AS'yi 
because IS.Sa is defined from 5"^ by subtracting the entanglement entropy for the ground state. For 
example, at t = we find 

ttMP vr M 

The time-evolution of ()3.7p is plotted in FigJSJ Notice that ASa respects the time-reversal 
symmetry. When I > a, the entanglement entropy initially grows with the time for t > and 



reaches the maximum at t = \/P — o? as in the left graph in Figl5j Later, it decreases as 

When I < a, the entanglement entropy always decreases for t > 0. At late time we have the 
behavior ()3.9p . Notice also that the width At of the peak around t ~ \/P — o? is estimated as 



At ~ a when q << L These results can be intuitively understood because at t = ^JX^ — the 
particle in AdS passes through the minimal surface ()3.5p . Refer to the Fig l2] again. In the dual 
CFT, they can be naturally understood if we remember the excitations propagate at the speed of 
light as we will discuss later. 



Notice also an interesting property that the height of the peak at t = — o? stays constant, 
given by Sa = iJ^Ri and thus does not decrease under the time-evolution. This can also be seen 
from the right graph in Figl5j This is rather different from the behavior of energy stress tensor 
studied in the previous section (see Figl3]). 

It is also intriguing to shift the center of the subsystem A relative to the trajectory of the falling 
particle in the (xi,?/2) plane. We take the minimal surface in the pure AdS to be 



P-{xi-0^-xl (3.10) 

where ^ is the distance between the center of A and the falling particle. 

We plotted the results of AS^ in Figl6] for different values of ^. In general one may notice 
that the peak is broadened so that it is spread over / — ^<t<Z + ^. In the gravity dual, this 
is easily explained from the propagation of gravitational waves from the falling particle. This also 
qualitatively agrees with the results of local quenches in two dimensional CFTs [S] . 



12 



Figure 5: The plots of ASa at ^ = for d = 3. The left and middle one describe Sa as a function 
oft. We choose {a, I) = (1,5) (left) and {a, I) = (1,0.5) (middle), respectively. The right 3d graph 
expresses Sa as a function of t and I for a = 1. The horizontal coordinate corresponds to /. We 
took the range -10 < f < 10 and < / < 10. We set R = AGn = M = 1. 






Figure 6: The plots of ASa for d = 3 as a function of time t. The left, middle and right graph 
correspond to (a, /, .^) = (1, 5, 0), (1, 5, 2) and (1, 5, 4), respectively. We set R = AGn = M = 1. 



It is straightforward to generalize the above results to other dimensions. In d = 2 and d = 4, 
the expression of ASa (|3.7p for d = 3 is replaced witlj^ 



d 



d 



ASa 
ASa 



2Mla + M{f 



o? — t^) arctan 



SGnIRo 



ttM 



SGnIRoi 



2al 



-, 3 ia^ — + t^) arctan ( -k I 



(3.11) 

+ 4a/ 



(q,2_;2+j2)2 ' 2al 

We can confirm that the behaviors of ASa in d = 2 and d = 4 are very similar to that in d = 3 
and thus we will not show them the explicitly. 
At late time t >> I, they are approximated by 



d 
d 



2 : ASa ^ 
4 : ASa - 



SGNRt^' 
5GNRt^ ' 



(3.12) 



This shows that ASa decays like 



t at late time. Notice that the perturbative calculation 



is always justified at late time as the back-reaction to ja for a finite / clearly gets suppressed. 



In d — 2 the subsystem is chosen to be an interval —l<x<l. 



13 



3.3 Small Subsystem Limit: An Analogue of the first law of thermodynamics 

In the small limit of the subsystem size /, we can trust the perturbative results of AS a ()3.7p and 
p.lip . This is because the surface ja is situated near the AdS boundary and therefore the deviation 
of the metric from the pure AdS is very small [29]. In asymptotically AdS backgrounds which are 
static and translation invariant, a relation which looks like the first law of thermodynamics has 
been found in the small size limit of A 1291: 



T,ffASA = AEA, (3.13) 

where Ea is the energy in the subsystem A given by Ea = dx'^~^Ttt- The effective temperature 
is defined by Te// = where the constant ca only depends on the shape of the subsystem A and 
is independent from the details of the CFT we consider. When A is a d — 1 dimensional ball with 
the radius / as we choose in this paper, we have 

T.„ = =^. (3.14) 

Thus it is intriguing to see if this relation (j3.13p holds in our time-dependent and inhomogeneous 
setup of local quenches. We can calculate the energy density from (j2.24p . and AS a from (|3.7p and 



()3.1ip in the limit / << yja^ + 1"^ as follows: 



d = 3: Ttt = —pz — ^^7-75— — rppr, AS a 



nGj^Rit^ + a^r ^"5G^i?(t2 + a2)4- (3-15) 
By using these expressions we can explicitly confirm the relation (|3.13p for any d. 

3.4 Large Subsystem Limit 

Before we go on, we would like to write down the result in the large subsystem limit / >> + a^. 
In this limit, we find from (13. 7p and (13. lip : 

d = 2: ASa = —pr'rnR, 
3t"^ 

d = 3 : ASa = —pr-mR, 

d = 4: ASa = ^^mR. (3.16) 

Note that these are time-independent. 
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4 Exact Holographic Entanglement Entropy for 2d Local Quenches 



We can actually find exact extremal surfaces (i.e. geodesies) in the AdSs case (d = 2). We take the 
subsystem A to be an interval at a constant time in the dual CFT2. In the holographic calculation, 
first we obtain a geodesic in the metric (j2.13p . which is asymptotically global AdSa. Then we map 
it into an asymptotically Poincare AdS metric by the coordinate transformation ()2.10p . 

Initially, we will assume M — E? < and thus the geometry ()2.13p has the deficit angle at r = 0. 
If we were allowed to take 9 to be 27t-^^=^ instead of 27r, the geometry gets smooth. Later we 
will come back to the case M — > i.e. the BTZ black hole. 

In the near AdS boundary limit, the map between the point (tqo, ^00 ^00) in (|2.13p and the 
point {t, Zoo, Xoo) in the Poincare AdS is given bj|f| 

2Rt 



tan Too 
tan 9nr, 



i?2e/3 + e-/3(x2„ -i2))' 

2it:xoo 



1 /„n o 1 



'R'^lo + i (e-f^i^lo - t') - R'e^f. (4.1) 

Notice that Zqo is interpreted as the UV cut off (or lattice spacing) in the dual CFT. In the above 
expressions, we chose the range of Too and 9oo to be [—tt,-k]. 
We can specify the geodesic 7a in ()2.13p by 

T = T{e), r = r{e). (4.2) 

Its length 1 7^ I reads 



|7A| = / d9^r^ + ^,^%_^ r'^ " (-^ + " M)r'\ (4.3) 



The minimal length condition is summarized as 



d9 r'^ + R'^ -M 
dr r 



JA^r'^ + (S2r2 - l)(r2 + i?2 _ m), (4.4) 
d9 R 

where A and B are integration constants. 
4.1 Symmetric Intervals 

Consider the case where the subsystem A is given by an interval — Z < x < Z at time t. The 
excitations of the local quench occur at x = because the falling particle is situated at x = 0. In 



^To fix the signs of Too and Soo, we need to go back to the original coordinate transformation (|2.10p . 
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the global AdS, this is mapped into an interval — ^oo ^ ^ ^ ^oo at a constant time slice r = Tqo, 
where we can assume < 0oo < without losing generality (see FiglT]). Since we can assume r is 
constant on 7^1, the extremal surface condition ()4.4p gets simplified to 

^ = -^y^{r^ + R^-M){r^-r^,), (4.5) 

where r = r* is an integration constant and r^, is the minimum value of r on 7^, i.e. the turning 
point. 

4.1.1 Case 1: M < R^ 

Let us first assume < Ooo < f- The curve with the minimum length is given by the geodesic 
which takes the angle range —6*00 < < Ooo (see the left picture of Figl?!). The point ^ = 
corresponds to the turning point r = r^, . Thus we finqj 

°° , Rn 
ar- 



ryf{r^ + R^ - M){r'^ 

R 



2^R^ - M 

It is useful to rewrite this into 



vr [R^-M-rl 

h arcsm — 7^ ^ 

R^-M + rl 



(4.6) 



\ R = rl + R^-M - ('-'^ 



Finally, the HEE is given by 



R 



Sa = -777^- / dr 



2Gn Jr, y^(r^ + R2_M){r^-r^) 

log /^o =^ ■ (4-8) 



2Gn v^i?2 _M + rl 

In the case ^oo > ^, on the other hand, the curve 7^ is given by the geodesic which takes the 
angular range ^oo < |^| < vr (see the right picture of Figl7|). Thus its turning point is at = vr. This 
guarantees the basic property of entanglement entropy written as Sa = Sb for pure states, where 
B is the complement of A. Thus for Q^o > §, we find the correct HEE from (|4.8p by replacing Q^o 
with vr - in (lO) . 

4.1.2 Case 2: M > R^ 

For M > R? , the geometry ()2.13p becomes the BTZ black hole. In this case, the basic result can be 
obtained from the previous one by the analytic continuation \/R'^ — M — )• i^J M — R^: the relation 



*In this paper, the function arcsin(a;) takes values between — |- and ^ . 
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0<6'., <nl2 



>nll 




Figure 7: The sketch of 7^ in the asymptoticahy global AdS background (|2.13p . The left and right 
picture corresponds to < Qoo < 7r/2 and 0oo ^ 7r/2, respectively. 

()4.7p is replaced with 

- i?2 + M 

= JL ^ f4 9) 

rl + R^-M' ^ ^ 

For < ^00 < f, the holographic entanglement entropy Sa can be found from the same expression 
(14. Sp . which is denoted by ^^(^oo)- For ^oo > f > "S'a is given by 5^1 (vr — ^oo) as in the M < i?^ case. 

In this calculation we implicitly assume that we are considering a star solution and that outside 
of the star is described by the BTZ black hole solution (j2.13p . Thus this is dual to a pure state in 
the CFT and indeed our construction satisfies Sa = Sb- 

4.1.3 Final Results 

The final result is plotted in FiglSl For a small M the exact result nicely agrees with the one from 
the perturbation theory (j3.1ip . Even for a large M, the agreement is very good except the peak 
at t ~ — q;2. Indeed, we can analytically show that in both the limit t — t- and t = 00, Sa 
approaches to 

5a(0) = SA{t = 00) = I log — , (4.10) 

where we employed the well-known relation c = [30] . This reproduces the well-known result 
of the entanglement entropy for ground states in CFTs [STJ [321 [2] . Remember that z^o corresponds 
to the UV cut off (lattice spacing) in the dual CFTs. 
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Figure 8: The exact plots of AS a as a function of time t(> 0). The left, middle and right graph 
correspond to M = 0.1,0.95,5, respectively. The blue and red curve correspond to the result 
from the exact analysis and perturbation theory, respectively. They coincide almost completely for 
M = 0.1. The dotted black one for M = 5 corresponds to the thermal local quench given by (j4.1ip . 
We set i? = a = 4Gn = 1 and / = 5. 

4.1.4 Thermal Local Quenches 

In the previous analysis for M > R^, we obtained Sa by replacing a black hole with a star in order 
to have a gravity dual of a pure state, which is usually assumed in the study of quantum quenches. 
If we deal with a falling black hole instead of a massive star, we will obtain Sa as follows 



of black hole entropy and arises because 'Ja wraps a half of horizon. The topological condition of 
J A of HEE requires that the subsystem A should be homologous to 7^. Therefore in the presence 
of the black hole horizon, we cannot simply change the geodesic from the one passing through ^ = 
to the one passing through 9 = tt. If we want to deform ja in this way, we need also to wrap ja 
on the horizon as its disconnected part. We plotted Sa for this thermal local quench in the third 
graph of FigJH 

4.2 General Formulation 

Now let us extend the previous analysis to more general setups where the subsystem A is given 
by an arbitrary interval l^^^ < x < l^"^^ at the time t. The surface ja is defined by the geodesic 
curve whose two end points are given by [x, t, z) = {l^^\t, Zoo) {i = 1, 2) in the Poincare AdS3 and 




(4.11) 



where min{a, b} denotes the smaller one among a and b. The term 



2Gn 



corresponds to a half 
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equally by {T,6,r) = (t!^ ,0^ ,r^oo) in the global AdS. Note that their relations are given by: 

tanrW = — 

°° ^2e/3 + e-/3((/«)2_i2)' 



tan 0(1) 



.(•0 



(4.12) 



Then the HEE reads 



Sa 



R 



AG 



N 



1=1 



JO 



dr 



Br 



^^2^2 + (^2^2 _ l)(^2 ^J12_M) 



(4.13) 



By integrating (j4.4p we find 



U(2) _^(1)| 
I ' oo ' oo I 

1^(2) _ ^(1)1 



R 



VR"^ - M 
R 

y/R'^-M 



IT 

— + arcsin 



vr 

— + arcsin 



B^{M - R^) + - I 



(S2(i?2 _M) + A^- if + 4^2(^2 _ M) ^ 

B^R^ -M) + A^-l 
{B^{R2 - M) + A2 - 1)2 + 4^2(^2 _ M) ^ 



(4.14) 



If we assume < — 6*1^^1 < vr, the HEE is computed as follows: 



Sa 



R 



log(r(^) •r(^)) + log 



AB' 



{B^{R^ - M) + A2 - 1)2 + 4^2(^2 _ M) ^ 



R 



2 cos 



log(r« •rg))+log 



R 



2 cos 



R'^ -M 



.(4.15) 



In the case A^oo = — o2d\ > vr, we need to replace A^oo — ?• 27r — AO^x, as we did in the previous 
subsection. 

So far we assumed M < R? . If M > i22^ then the solution (j2.13p is a BTZ black hole without 
any deficit angle. However, our analytical calculation done before still holds via the analytical 
continuation as in the previous case. 



4.3 An Interval with An Excited End Point 

Especially, let us focus on the case = (0,/). In this case, one of the end points of A i.e. 

X = is excited by the local quench. The result is plotted in FigJH It is easy to see that ASa is 
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Figure 9: The exact plot of AS a as a function of t for Z*-^^) = (0, 5). The left, middle and right 
graph correspond to M = 0.01, 5, 500, respectively. We set R = a = AGn = 1- 



non-trivial during —l^t<l,as expected from the causality argument. This qualitatively agrees 
with the CFT result in [8]. As M gets larger, the HEE looks like a step function. 
In the late time limit t — )■ oo, Sa approaches to the result of the ground state 



SA{t = oo) = - log- 

•J •2oo 

On the other hand, at the time i = 0, we find the following result for I » a 



SA{t = 0) = SA{t = + 6 l°g I ^2 _ ^ 



Notice that we always have SAit = 0) > ^^(t = oo) 



• sm 



WR'^ - M 
2R 



(4.16) 



(4.17) 



4.3.1 Semi Infinite Limit I — t- oo 

It is useful to consider the limit I » t » a for = (0,/). We find 

-(2) 



2ta 



(1) 

' OO " ^5 ' OO ' 

^(1) = 0(2) ^ ^ _ ^ 



^ (1) 

OO — 



Rt^ 

2aZo 



.(2) 



RP 

2az„ 



Finally the HEE is found to be 



^ c , Rtl c , 

Sa-i: log h - log 

6 Zazoo o 



R"^ -M 



VR^ - M , , A . 
sm I (jx — ajt) ■ sm 



R 



C - / C - t c , 
~ - log h - log h - log 

2^oo 



R 



3" ° Zoo ■ 6 " °a 6 y - M 
In this way, we obtain the behavior 



sm 



'K^JR^ - M 
R 



S 



AdS 



c , t c . 
~ ^ log - + - lo, 



(4.18) 



ay/R^ - M 

m 



(4.19) 



(4.20) 
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Figure 10: The approximation of ASa by the function logt+const. for = (0, 1000). The 

thick black curve is the plot of the exact result, while the dashed one is its approximation by 
logt + 0.1. We choose M = 0.5 and set R = a = 4:Gn = 1- They nicely agree with each other 
when Q << t « I. 

for the time evolution. Refer to FigJlOlfor the numerical confirmation. 

On the other hand, the local quench induced by joining two half lines leads to the following 
behavior [8] 



where e is the cutoff for the process of local quench introduced in [8] and is analogous to a in our 
model. 

Note that both have the common important property that Sa increases logarithmically after 
local quenches. This is a characteristic feature of local quenches in two dimension. To be more 
precise, the coefficient of logt is different between (|4.20p and (14.2ip . Indeed, it is known that this 
coefficient depends on how we locally excite the system as found in the example analyzed in [7]. 
The difference of the coefficient in front of log / is easy to understands because in the case (j4.20p 
the system was already joined before the quench, while in the other case (14.2ip . the system was 
originally disconnected. We will give a more detailed explanation of the two difference behaviors 
(I4.20j) and (I4.2ip in section 6.2 using the tensor network description. 

4.4 HEE for General Intervals 

Now we would like to turn to the HEE Sa for a subsystem A defined by a general interval 
Equally we can think that A is parameterized by its width I and the position of its center ^ via 



We can explicitly calculate <S'yi(/,^, t) by using the formula (I4.15p . We plotted ASa as a function 
of ^ for fixed values of / and t in Figllll 

First, we notice that for a small / (the left graph of Fig llip . the initial peak at ^ = divides 
into two peaks under time evolution, keeping total sum of the heights conserved. This result looks 




(4.21) 



/(i)=^-//2, l'^'^^ =^ + l/2. 



(4.22) 
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Figure 11: The plot of AS a as a function of ^. We set I = 0.1 (left), I = 3 (middle) and I = 10 
(right). In each graph, the blue, red and black curve describes AS'yi(/,^,t) for t = 0, t = 5 and 
t = 10, respectively. We choose M = 0.9 and set i? = q = 4GAr = 1. 

very similar to the energy density Tu in two dimensional CFTs (see the left graph in Figj3]), where 
the energy conservation holds. Indeed, we can easily understand this coincidence from the first law 
relation (|3.13p in the small / limit, which can be rewritten as 

ASA{l,^,t)c^^Tu{C,t). (4.23) 
Therefore we can conclude that in the small I limit, the integral 

/oo 
d^ASA{l,tt), (4.24) 
-oo 

does not depend on the time t. 

Moreover, for / > a, we find that AS a at t = has a peak at |.^| = 1/2. In the gravity dual, 
this is easy to understand because the geodesic 7^ comes very close to the massive particle as 
explained in Figl2j At a later time t > 1/2, the hump propagates at the speed of light, conserving 
its form. This occurs because "ja is now far from the falling particle and its shock waves propagates 
without changing its form. Therefore we expect that the integral (14.24p will approach to a constant 
value at late time t > 1/2. Indeed, we can confirm this numerically as shown in the FigJ12l This 
behavior roughly tells us that the long range entanglement will be generated at late time, while the 
short range one will be not. We will study more carefully such a dynamical structure of quantum 
entanglement in the next section, which will provide clear explanations for all of the results in this 
section. 



5 Entanglement Density 

In this section, we will introduce a new quantity which we call entanglement density. It can extract 
an essential structure of quantum entanglement from Sa in two dimensional quantum field theories. 
We will apply this to our holographic local quenches in AdS3/CFT2. 
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Figure 12: The numerical plot of / (i^AS'^(/, ^, t) as a function of time t > 0. The blue, red 
and black graph correspond to / = 1,5, 10, respectively. In the small / limit, this quantity should 
approach — 0.262 as follows from (|4.23p . We choose M = 0.5 and set = q = ^Gn = 1. 

5.1 Definition of Entanglement Density 

Consider a general two dimensional field theory in a certain pure state and suppose that we can 
calculate its entanglement entropy Sa- The subsystem A is defined by the interval l^^^ < x < l^"^^ (or 
more simply [Z^^^/^^^]), where x is the space coordinate of the field theory. Just for convenience, 
we will treat Z^^^ and /^^^ to be independent parameters for a while. We would like to ask how 
to extract a structure of quantum entanglement explicitly from the entanglement entropy Sa for 
arbitrary intervals. 

For this purpose, we take into account only the two body entanglement and define the entan- 
glement density n{l, ^, t) such that this counts the number of entangling pairs (or bits) between the 
two points X = ^ - l/2{= /(^)) and x = ^ + l/2{= P^). Ima gine that we discretize the field theory 
into a lattice theory such as spin chains. In each point, we assume that there are several spins. 
Then n{l,^,t) counts the number of entangled pairs of spins which are located at x = — 1/2 and 
x = ^ + 1/2. Therefore the parameter / describes the range of the entanglement and ^ does the 
position of the center of the entangled pair. 

Now we approximate Sa by summing all entangling pair as follows (see Fig llSp : 



Sa = dx 
//(I) 



OO POO 

dw n{w,x — w/2,t) + / dw n{w,x + w/2,t) 
-iw Jm-x 



(5.1) 



By taking a derivative with respect to Z*-^-* with l^^^ fixed, we obtain 

dSA _ 
dm 



dw n [ w 



li^^-lt)+ rdwn{w,li^^ + ^,t] 



-£ydxn(^/(2)-x,| + ^,t^. (5.2) 

Finally by taking the derivative by we find 

a^S, lO^S, O^S, ,„(,<^)-«>..^^^l±ii!l,(U2„(«,*,, (5,3) 
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n(w,x-w/2,t) 
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n(w,x+w/2,t) 



Figure 13: A sketch of expression of Sa as a sum of entanglement density. Hooks with the arrows 
on both end points describe the entangling pairs between two points. 

In this way, we can extract n{l, ^, t) if we know Sa as a function of Z^^^ and t. 

For a ground state of a CFT2 with the central charge c, Sa is given by (|4.16p and thus the 
entanglement density is found to be 

ncFT{l,tt) = -^. (5.4) 

Notice that the entanglement density is free from the UV divergences. Thus, when we talk about 
a deformation of a CFT2 as in the present example, it is sometimes useful to normalize the entan- 
glement density by multiplying l^. 

5.2 Strong Subadditivity and Positivity of Entanglement Density 

Consider a two dimensional field theory in a generically excited state. We introduce three subsys- 
tems Ai,A2 and A3 so that they are defined by the intervals [/i,/2]; [^2)^3] and [Z3,/4], respectively, 
such that h < I2 < I3 < h- 

The strong subadditivity [33] (refer to [34] for its application to field theories and to [S^ for its 
holographic proof) is given by the inequality 

Sa + Sa 2UA3 >Sa + Sa,. (5.5) 

In our setup, this is equivalent to 

S{h,k) + Sih, k) > S{l2,h) + S{h,h), (5.6) 

where S{p,q) is the entanglement entropy Sa for the interval [p,q]. 
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In particular, we can take the limit 61^2 ~^ +0 with I2 = h + ^1 and I4 = 13 + 62- Then, the 
strong subadditivity (I5.5p is rewritten into 

This proves that the entanglement entropy in (I5.3P is positive. 

5.3 Conservation Law 

It is useful to study how the total integration of the entanglement density behaves. To make the 
analysis simple, we assume that the x coordinate can be compactified on a circle with the length L. 
The infinitely extended system we are mainly studying in this paper can be obtained in the limit 
L — )■ 00. Since we have the obvious relation n{l,^,t) = n{L — l,L — ^,t) due to the periodicity, 
we can restrict the integration of / to < / < L/2. We will focus on the increased amount of 
entanglement entropy and entanglement density compared with the ground state, denoted by I^Sa 
and An, respectively. 
Thus we find 

j-L /■L/2 j-L i-lW+L/2 

dU dlAn{l,^,t)= dl^^n dl^^^ An{l,C,t) 

Jo Jo Jo Jiw 



2 Jo JiW 5/(1)5/(2) 



L 







d 



d 



ASa(x,/(i)+L/2) A5^(x,/(i)) 

x=lW OX 



x=lW 



(5.8) 



dx 

We can show that the first term in the third line of (j5.8p vanishes. This is because by employing 
the identity Sa = Sb for pure states and the periodicity of the circle, we can easily show 

r dy [A5a(2/ + V2, y + L/2) - ASAiv - S/2, y + L/2)] = 0, (5.9) 
Jo 

for any S. Moreover, the other term in the third line of (j5.8p vanishes because the first law (j4.23p 
tells us 

ASA{y + 6/2, y) - ASA{y - 6/2, y) ~ |r« • 6^ (5.10) 

in the 6^0 limit. Remember that our first law assumes that there is a UV fixed point in the field 
theory we consider. 

In this way, we proved that the total integration of entanglement density is constant. In other 
words, the total number of entangled pairs is conserved: 

rL fL/2 



[ dC [ dl Anil,C,t) = 0. (5.11) 
Jo Jo 
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In general, the negative contribution comes from the UV region where I is very small. Indeed, 
by using the first law ()4.23p . we can generally show the following result in the limit / — >• 0: 

limAn(/,e,t) = -^r«(C,i). (5.12) 

5.4 Entanglement Density for Global Quenches 

To understand better the meaning of entanglement density, we would like to consider the global 
quenches in two dimensional CFTs with central charge c before we examine the local quenches. As 
shown in [6], the increased amount of entanglement entropy for the subsystem with the width / 
behaves like 

ASa = c - Am-t {e<t<l/2), ASa = ^ - Am-l (i > //2), (5.13) 

where Am is the energy gap which creates the global quench times a numerical constant and is 
proportional to the effective temperature at late time; e is cut off scaled of the quench process, 
which is of order 6m~^. 

In this case, the entanglement density (|5.3p reads 

c ' Aui 

An{l, t) = ■ S{t - 1/2) + An{l, t)uv, (5.14) 

where the last term denotes the UV contribution which is non-trivial only for / of order e ~ Am^^ 
or smaller. This UV term is negative such that the conservation law (j5.1ip holds. The delta 
functional term in (|5.14p shows that the entangled pairs have the range 2t at time t. We sketched 
this behavior in Fig|14l This behavior can be understood that the entangled pair moves in the 
opposite way at the speed of light, being consistent with the interpretation in [6]. 

In this way, for the excited states produced by quantum quenches, An{l,^,t) tends to positive 
for large /, while it becomes negative for small / due to the conservation law and the first law as in 
(|5.12p . Therefore we learn that an excitation of a quantum system corresponds to breaking shorter 
range entanglement and creating longer range entanglement. 

5.5 Entanglement Density for Holographic Local Quenches 

Now let us come back to the analysis of holographic local quenches. We consider the difference from 

the ground state and would like to calculate An{l, ^, t) from A5yi(/, t) by using the formula (j5.3p . 

The results are plotted in Figll5[ We can confirm that the total (properly normalized) density 

/^n(/,^,t) = 1 + /^An(/,^,t) is always positive as expected (we set c = 6 in the plot). 

More precisely this is the length scale beyond which we can approximate the excited state by the boundary 
state, which preserves the boundary conformal invariance. See [14] for more details. 
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Figure 14: A sketch of the entanglement density as a function of I under a global quench. 

First notice that the positive peaks in FigllSI are all delta-functional. This is because our 
gravity backgrounds, where we calculated the HEE, are singular due the presence of the deficit 
angle, corresponding to the point particle limit. This delta-functional peak appears at 

e = ±^/P/4-t^ -a'^, (5.15) 

which is when the massive particle is on top of -ja (see Figj2]). However, it is easy to imagine a 
regularization by replacing the point particle with a finite size object, which will replace a delta- 
functional peak with a smooth peaks with a finite width (called ij). 

Employing this regularizationlf] the entanglement density n(/, ^, = n{P^-l^^\ (/(2) + ;{i))/2,t) 
is plotted as a function of l^^^ and /^^^ in FigJ16l for t = and t = 2. It is clear from this graph 
that the peaks exist along the curve (f05D or equally /(i) • = 

— — t'^ as expected from the 

condition that the massive particle passes through 7^. Moreover, An takes the largest value when 

1(2) ^ _^{1) ^ ^/^2^, 

These behaviors can be systematically understood by considering the time-evolution of each 
entangled paiiQ. At time t, the pair of the points x = ^ — \/^^ + + and a; = ^ + \/^^ + + 
are entangled for any values of ^. Especially the pair at ^ = (i.e. the one between x = —\Jo? + 
and X = \/a^ + i^) possesses the strongest entanglement, which we call the dominant entangled 
pair. More generally the entanglement is enough strong for the entangled pairs with |,^| < a. 
Therefore they follow the time evolutions in an obvious way, as sketched in FigJ171 At earlier time 
®More explicitly, we replaced the derivative g^fi)g^(2) with the second order difference Sa{1^^^ + rj/2,l^^^ + v) ^ 

^Another holographic realization of entangled pairs refer to [36] based on the AdS/BCFT proposal [37] . 
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t << ^ the pair moves slowly. At late time t >> ^, they move at the speed of light in the opposite 
direction. This time evolution is also sketched with the energy density in FigHl 

We can explain the qualitative behaviors of HEE calculated in the previous section by referring 
to this structure of quantum entanglement in our holographic local quenches, based on the entangled 
pairs. For example, let us explain the result described in Figll2[ First assume t < 1/2. Then if 
we sweep a length / interval A by shifting its center ^, Sa_ cannot detect the entanglement pairs 
generated by the local quench when t — //2 < ^ < — t + //2 because the dominant entangled pair is 
completely included in A for these values of S,. This missing contribution is proportional to 1/2 — t 
and therefore Sa increases linearly until t = 1/2. At late time t > 1/2, for any ^, the interval 
A always include only one of the pair and therefore it fully contributes to Sa- Thus the integral 
J (I^ASa no longer depends on the time t. These explain the behavior found in FigJ12i 

One may think that any long range entanglement should not exist at t = because the massive 
particle is situated at deep UV region z = a. However, this speculation is based on a too naive 
UV/IR relation of AdS/CFT. Indeed, if we consider a two point function (O(x)O(y)) for a certain 
operator with a large conformal dimension, then we can evaluate this from the geodesic distance 
in AdS. The same geodesic distance appears in the calculation of HEE and therefore we find that 
a non-trivial result for two point functions can obtained when x = 0, where the geodesic 7^ passes 
very close to the massive particle as in right picture of Figl2j Thus this consideration of two point 
functions also support the long range entanglement even at t = 0. Another way to see this is to 
note that the energy density Tu computed in (I2.24p decays only slowly as a;~^ in the long distance 
limit X —7- 00 for any non-zero a. 

In this way, we learned that the long range entanglement is non-trivial even at t = in our holo- 
graphic model. In ideal models of local quenches we may not expect any long range entanglement 
just after the quench. However this is not a major problem because as we mentioned the dominant 
entangled pair is the one between x = —\Jo? -|- and x = Va^ + 1^, which is short range (~ a) at 
t = 0. 

In summary, we can understand the time evolutions of entanglement entropy 5^ in terms 
of those of entanglement density n. In our holographic local quench, the dynamical behavior 
of n is described by the evolution of the entangled pairs as in Figll7l At late time, the range 
of entanglement increases at the speed of light. As the time evolves, short range entanglement 
disappears and long range one is generated so that the total number of entangled pairs is conserved 
as in (jS.lip . This process looks a sort of the decoherence phenomena. 
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Figure 15: The numerical plots of PA7i{l,^,t) as a function of the center position ^ at t = 
(left graph) and t = 2 (right graph). The blue, red and black graph correspond to / = 1,2,10, 
respectively. We choose M = 0.5 and set R = a = ^Gn = 1- Notice the relation Pn{l,(^,t) = 





Figure 16: The numerical plots of An{l, ^, t) as a function of /^-^^ = ^ — //2 < (horizontal axis) and 
/{2) = ^ + //2 > (inward axis). We choose the regularization parameter as rj = 0.4. The left and 
right plot corresponds to t = and t = 2, respectively. We only plot those points where An(Z,^,t) 
is positive. We choose M = 0.5 and set i? = q = AGn = 1. 

6 Quantum Information, Thermodynamics and Gravity 

In this section we employ the results found in the previous sections to develop some insights on 
connections between the quantum information and gravity via the AdS/CFT correspondence. 

6.1 Total Amount of Quantum Information 

In the CFT side, the entanglement entropy Sa is interpreted as the entropy for an observer in A, 
assuming that the observer is not accessible to the region B. In order words this measures the 
amount of quantum information inside B for such an observer. Therefore ASa can be regarded as 
the amount of quantum information in B which is induced by the local quench, such as the data 
of excited particles etc. 

Notice that the standard thermal entropy for a mixed state, which is dual to a black hole entropy 
via AdS/CFT, corresponds to the information which is not accessible to the observer by any means 
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Figure 17: A sketch of time evolution of entangling pairs. The blue thick curve describes the 
dominant entangled pair ^ = 0. Around this, other entangled pairs are gathered with the width of 
order a. 

(except for making microscopic observations of its heat bath). On the other hand, the entanglement 
entropy Sa measures the amount of information which is possessed by B and is accessible to the 
observe in A after all possible experiments are done. 

Therefore, the holographic entanglement entropy allows us to calculate such a novel quantity, 
i.e. the amount of information ASa, for any massive object in AdS, which is dual to a certain 
excited state in the CFT. Apply this idea to our falling particle in AdS. The standard principle in 
AdS/CFT tells us that the size / of the localized excitations induced by the local quench is estimated 
as I ~ \/t^ + a^, which is the value of the coordinate z of the falling particle. Therefore it is natural 



that we choose the subsystem A (or equalljo B) to be the ball with the radius + a^. Indeed, 
we can see from our perturbative results (j3.7p and (jS.lip that ASa takes the maximal value when 



I = yjt'^ ^ o? as a function of I and takes a time-independent constant value. Therefore we call it 
^^max_ 'pj^jg ig evaluated as the following simple form for any dimension d: 

^^max = Cs-mR = Cs- A, (6.1) 

where Cs is an order one constant. Explicitly, in each dimension d, we find = 2 (d = 2), 
Cs = it/2 {d = 3) and Cs = 4 [d = A). We used the well-known relation between the mass m of a 
particle and the conformal dimension A of its dual operator [10] assuming A >> 1. For example, 



Notice that since we are considering a pure state we always have Sa = Sb and therefore the distinction between 
A and B is not important. 
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in the AdSa case the exact calculation (|4.8p leads to at / = + a^: 



c 



A5- = ::iog I sin MEIE] . (6.2) 

This shows that AS""°^ is a constant which is independent from I. Assuming the range 1 << A << 
c we fincjf 

^gmax ^ 2A. (6.3) 

In this way we find the number of miscrostates of the massive particle is estimated as follows 
in any d: 

# Microstates ~ e^"^"""" = e'^"'^. (6.4) 

This is qualitatively consistent with the well-known Hagedron growth of degeneracy of the states 
with a large conformal dimension A > > 1 \25\ [38] . Remember that this relation ()6.4p holds only 
when the backreaction of the massive particle is small. For example, this requires A « c in d = 2 
as can be seen from (j6.2p . 

If we remember that the energy E of localized excitations made of radiations in the CFT is 
given by ()2.5p . then we find the simple relation between the amount of information of this excited 
lump and its the total energy: 

^gmax ^E-a. (6.5) 

Remember that a denotes the size of the lump at the initial time t = 0. This provides us a 
prediction that in the large strong coupled gauge theories, the amount of information possessed 
by such a radiation 'fire ball' is given by its energy E times its linear size a. It would be very 
intriguing future problem to extend this estimation to other objects like the fundamental strings 
or D-branes in AdS/CFT in order to see how much the relation (jO.Sp is universal. 

Finally, we would also like to comment on a possible interpretation of AS"""^ from the gravity 
side. Since the Unruh temperature of the falling massive particle reads Tu = at t = 0, we can 
rewrite this relation (16. 5p into 

Tu ■ AS'"''^ ~ E, (6.6) 

which looks like the thermodynamical first law. We can speculate that this is the thermodynam- 
ical relation which is understood as that for a Rindler observer, whose acceleration is a~^. This 
interpretation makes sense near t = where the energy from the viewpoint of the Rindler observer 
coincides with E. 



^If we assume a very large dim. operator A >> c in the local thermal quench, then it is dual to a BTZ black 
hole. In this case we find from (|6.2|l AS"""" ~ = |\/3cA. This is a half of the black hole entropy and agrees 

with the Cardy formula. 
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6.2 Entanglement Renormalization and Origin of Gravitational Force 

A helpful framework to study the connection between quantum entanglement of a given quantum 
many-body system and its holographic geometry is the entanglement renormalization (or called 
multi-scale entanglement renormalization ansatz, MERA) [39j. Indeed, it has been conjectured 
in |40] that the geometry described by a tensor network of a quantum critical system in MERA 
describes the AdS spacetime dual to this critical point. In [51], this correspondence has been 
studied in the continuum limit (i.e. field theory limit), called cMERA [32] and a candidate of 
holographic metric in terms of purely field theoretic data has been proposed. Refer also to [43| for 
recent developments in this subject. 

Motivated by this it is intriguing to see how we can describe the time evolution of local quenches 
from the viewpoint of the entanglement renormalization. A general framework to describe a quan- 
tum state graphically is called tensor networks (see e.g. reviews [44]) and has been recently used 
to prepare a candidate of approximate ground states based on the variational principle, optimized 
by a number of parameters. 

The tensor network for a CFT is called MERA and the network for its ground state is described 
in FigJlSI It describes the coarse-graining procedure of a given quantum system such as spin chains 
as we goes from the bottom (UV) to the top (IR), by combining two sites (or spins) into one, which 
looks like a tree in FigllSI To realize the entanglement structure of CFTs, we need to add so called 
disentanglers which are unitary transformations between two spins and which are written as the 
horizontal ladders in FigJlSI 

Let us define the number of disentanglers such that each of them carries the entanglement 
entropy log 2. Then we find that the number of disentanglers in each bond of the MERA network 
is given by |. Indeed, this reproduces the correct entanglement entropy for ground states in CFTs 
given by the well-known formula (j4.10p as explained in FigllSI Notice that we can estimate the 
entanglement entropy from the MERA diagram by encircling the subsystem A and counting the 
number of intersecting bonds (disentanglers). Though there are many choices of encircling curve, 
we can optimize the result by choosing the one with minimal number of intersecting bonds. This 
minimizing procedure nicely matches with the HEE formula (j3.1|) as pointed out in j40j . 

By comparing this estimation of the entanglement entropy in CFTs with (|5.4p , we can speculate 
that the number of disentanglers in a bond is roughly approximated by P ■ n{l, t) in more general 
tensor networks which do not correspond to ground states of CFTs. Here the averaging symbol 
means an average like ~ ^ Jq dl f^^il2 d^- The conservation law (jS.lip tells us that the total number 
of disentangles is also conserved. Moreover, we can roughly speculate the relation between the 
entanglement density and the holographic metric in the extra direction z: \/gzz{z, x, t) oc I -nil, ^, t) 
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Figure 18: The tensor network of MERA for a ground state of a CFT. The coarse-graining procedure 
is done from the bottom (UV) to the top (IR). The horizontal ladders represent the disentanglers, 
while the tree structure does the coarse-graining procedure, combining two sites into one. We also 
show the estimation of entanglement entropy when the subsystem A is an interval with length / by 
counting the number of bonds. 

with the identification / = z and ^ = x. It will be an interesting future problem to work out these 
relations more precisely. 

Now we would like to consider how to describe local quenches by using tensor networks. We 
consider both our holographic local quench and the original local quench in [8]. We argue that 
they are described by the left and right pictures of Fig ll9l respectively. In our model, there are 
disentanglers even in IR regions because we start with the ground state of a CFT on an infinite line. 
On the other hand, in the local quench induced by the joining two semi-infinite lines, there should 
be no entanglement between the two in the IR region. As in (|4.10p and ()4.16p . we can confirm 
that the behavior of Sa when t >> / is the same as that for the ground state. Thus the MERA 
structures in the UV region are also the same as that for the ground state in both cases. Finally, 
the massive falling particle corresponds to the insertion of many disentanglers. This description is 
consistent with the behaviors ()4.20p and (j4.2ip at late time. 

Finally we would like to ask what is the holographic origin of gravitational force. As we have 
seen, the propagation of entangled pairs is dual to the falling particle in AdS. This tells us that the 
evolution of the range of entanglement is dual to the gravitational force in AdS space. Therefore 
we can argue that the gravitational force is dual to a kind of decoherence which breaks the short 
range entanglement and tries to expand the range of entanglement. In MERA, this is described by 
the motion of the dense lump of disentanglers toward upper direction as depicted in FigJTOl 
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Figure 19: The tensor networks of MERA after local quenches in a CFT. The left picture corre- 
sponds to the local quench triggered by joining two CFTs on semi-infinite lines. The right one 
describes our holographic local quench induced by local excitations in the CFT. 

7 Conclusions and Discussions 

In the first half of this paper, we proposed a holographic model of local quench and study its property 
especially focusing on the time evolution of entanglement entropy. In our model, the local quench 
is simply described by a free falling massive particle in an AdS space. We performed a perturbative 
analysis of its back-reacted geometry and computed the holographic entanglement entropy (HEE) 
in any dimension as well as the holographic energy stress tensor. Moreover, we presented exact 
calculations of HEE for local quenches in two dimensions. We leave exact calculations in higher 
dimensional cases for a future problem. Since in our model we considered a massive particle without 
any charge, it will be another interesting future problem to extend our results to charged massive 
particles. 

In most of earlier works on local quenches, local excitations are generated by joining two semi- 
infinite lines (or spin chains). On the other hands, our holographic model corresponds to the setup 
where local excitations are induced in a CFT on an infinite line and does not involve any joining 
procedure. Therefore we found that details of the time-evolution of entanglement entropy are 
quantitatively different between them, while we find that the qualitative behavior agrees with each 
other. For example, the logarithmic time evolution ~ logt is common to both of them for a large 
subsystem. However, its coefficients differ by the factor two between them and we give an simple 
explanation from the viewpoint of entanglement renormalization. This is because there exists long 
range entanglement even just after the quench. Even though we were able to reproduce the basic 
properties of local quenches in our model, it is an intriguing future problem to construct a clean 
gravity dual without any long range entanglement. 
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In order to study the time evolution of quantum entanglement more clearly, we introduced a 
new quantity called the entanglement density for two dimensional field theories. This measures 
the density of entangled pairs between given two points. We showed that the strong subadditivity 
guarantees that this quantity is positive as should be so if we want to understand it as a physical 
density. Moreover, we found a simple relation between the entanglement density and the energy 
density in the small subsystem limit. 

It will be interesting to generalize this quantity so that we can incorporate the effect of three- 
body or higher entanglement. At least it is already clear that A^-body entanglement is related to a 
iV-th derivative of the entanglement entropy for a suitable subsystem. At the same time, it will be 
an important problem to extend the entanglement density into higher dimensional counterparts and 
study various holographic models. A special example of higher dimensional entanglement density 
is obtained by replacing the interval subsystems with strip ones. 

In the final part of this paper, we calculated the amount of quantum information carried by 
a massive object (localized excitations) in the dual CFT. We find a rather simple rule that this 
amount of information is given by the energy E of the object times it linear size a. This is expected 
to be a prediction which is true for large N strongly coupled gauge theories. It will very intriguing 
to see how this behavior is universal in more general cases. 
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A Perturbative Construction of Gravity Duals 

In this appendix, we will show a calculation of the back-reaction for a falling massive particle in 
AdS4 by solving ([221) and i^B) directly. 
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A.l Perturbation Theory 

We focus on d = 3 i.e. the AdS4 case and define x = xi and y = X2- We will set i? = 1 just for 
simplicity. We perform the Fourier transformation with respect to x, y and t. Using the symmetry 
X —7- —x and y — t- —y, we can write down the following ansatz of the perturbed metric: 

Sg.. = {klfiz) + g{z)) e^ik^-^+f^yy-'l 6g,y = k,kyfiz)e'^'^--+'yy"^'\ 

6gyy = {klfiz) + g{z)) e^i'^.-^+'^yy-t) . (A.l) 

The energy stress tensor for the falling particle described by z{t) = \fW^rO? after the Fourier 
transformation 



T^-(k^,ky,oo,z) = ^ j dxdydt e-'^'^^-^+'^yy-'^'^ ■T^''{x,y,t,z). (A.2) 



is given by 



^V^2 - -2 

a 



z^ 



T^* = 2im ■ — • sinfwv — q;2^ 
a 

^7 

r" = 2m • cos(wVz2 - a2) 



(A.3) 



and other components are vanishing. We employed the identity 

^/z^ — (T^ V ' 



\f z^ — a 

where we assumed z > 0. We also defined 'normalized mass' 

STrGArm 

It is straightforward to check that they satisfy the conservation law V^T'"^ = 



^^o^^^^ (A.5) 



A.2 General Solutions to Einstein Equation 

By evaluating the left-hand side of (12. 7p . in the end we find that all the solutions can be constructed 
as follows: First we define m(z) such that 

, \ A,. N , \ U'{Z) U"(Z) 

w{z) = A{1 -u:) + —u{z) - ^ + 

' u: 2 zu^ 2w2 ' 

_ v!{z)_ ^ ^u'{z) _ u"{z) _ ^u"{z) ^ 2u^^\z) _ -uW(z) 

z z^ct;2 2 z2a;2 za;2 2ct;2 ' 

, , , ^ 2n'(z) Au\z) , 4u"{z) 2n(3) u^^) 

Hz) = A. - -uiz) + ^ + ^ - uiz) - ^ + ^ - (A.6) 
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Moreover, we define y{z) by 

= (A.7) 

Then we can show that ()A.6P is a solution to (j2.7p if the fohowing equation is satisfied 

y^^\z) - {kl + kl- uj'^)y(z) - ^ [sin (u^z^ - a^) - ujVz^ - cos (uy^z^ - ^2)] =0. (A.8) 



It is clear that there are five integration constants. They are dual to the background metric 
perturbation (four parameters) and the boundary energy stress tensor (one parameter). 

A. 3 Special Solutions: a — )■ 

We can solve (jA.Sh in the particular case a — )■ where z(t) ~ \t\. We impose the boundary 
conditions such that there are no non-normalize deformations (i.e. the background metric is not 
perturbed) and that the metric perturbation does not blow up in the IR limit z ^ 00. This uniquely 
fixed the form of u{z) as follows 

""^^^ = ~ aa;4(fcg + fcg - (^^' + ~ + + ~ 

+ ^ ( (-4A:2 _ 4A;2 - 2lo^ + {kl + ^2)^2^2^ cos(2a;) - 2zLo{2kl + 2k^ + uj^) sm{zuj)) 

{kx + ''^yj 

+ (^2+^2)2 • 2e-V^^i^^a;6 (l + z^ k^^ + k^ - . (A.9) 

We also have the condition ^ = 0. These completely fix f{z),g{z),h{z) and w{z). 
Near the AdS boundary z — )• they behave 

, , 4r7ia;2 
i(j(z) ~ • z , 

M^l + kD^kl + kl-uj^ 

Am{kl + kl- 2u?) 3 

/(z) ~ , • z\ 

Mkl + kl)\kl + kl-u:'^ 



Am{kl + kl- ^2) 
g{z) ~ , • z-", 



'iaikl + kD^kl + kl-u:'^ 

/i(z) ~ — , • z^. (A.IO) 

3JA:2 + fc2-a;2 
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We can read off the tensor tab defined by (|2.17|) as follows: 

4m 

ttt 



ttx 



tty 



ArhkrUJ 



AmkyUJ 



txy 



4m{k^ + klkl - k^u? + aj'^k^) 
2,a{kl + klY ^ kl + kl - u:"^ 
Arhkxky{k1 + ky - 2a;^) 



3a{k'^ + k'^y^/kx + ky - co 



,2 



Am(k^ + klk"?, - kluj"^ + uj'^k'}) 
tyy = — ^ ^ (A.ll) 

Mkl + kl)\lkl + kl-u:^ 



y, Y , 

The physical holographic stress tensor Tab is calculated from (I2.17P via the formula (I2.18p . In 
the end we will find that (jA.lip reproduces the result (I2.20p at a = 0. For example, if we take the 
Q — )■ limit of ()2.20p . we can show 

by using 

lim ^ = ^6{x). (A.13) 

By performing the Fourier transformation, we find 



This precisely agrees with (jA.lip . 
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